The endosymbiont of Epithemia clementina is specialized for nitrogen fixation within a photosynthetic eukaryote.

Epithemia spp. diatoms contain obligate, nitrogen-fixing endosymbionts, or diazoplasts, derived from cyanobacteria. These algae are a rare example of photosynthetic eukaryotes that have successfully coupled oxygenic photosynthesis with oxygen-sensitive nitrogenase activity. Here, we report a newly-isolated species, E. clementina, as a model to investigate endosymbiotic acquisition of nitrogen fixation. We demonstrate that the diazoplast, which has lost photosynthesis, provides fixed nitrogen to the diatom host in exchange for fixed carbon. To identify the metabolic changes associated with this endosymbiotic specialization, we compared the Epithemia diazoplast with its close, free-living cyanobacterial relative, Crocosphaera subtropica. Unlike C. subtropica, in which nitrogenase activity is temporally separated from photosynthesis, we show that nitrogenase activity in the diazoplast is continuous through the day (concurrent with host photosynthesis) and night. Host and diazoplast metabolism are tightly coupled to support nitrogenase activity: Inhibition of photosynthesis abolishes daytime nitrogenase activity, while nighttime nitrogenase activity no longer requires cyanobacterial glycogen storage pathways. Instead, import of host-derived carbohydrates supports nitrogenase activity throughout the day-night cycle. Carbohydrate metabolism is streamlined in the diazoplast compared to C. subtropica with retention of the oxidative pentose phosphate pathway and oxidative phosphorylation. Similar to heterocysts, these pathways may be optimized to support nitrogenase activity, providing reducing equivalents and ATP and consuming oxygen. Our results demonstrate that the diazoplast is specialized for endosymbiotic nitrogen fixation. Altogether, we establish a new model for studying endosymbiosis, perform a functional characterization of this diazotroph endosymbiosis, and identify metabolic adaptations for endosymbiotic acquisition of a critical biological function.

The endosymbiont of Epithemia clementina is specialized for nitrogen fixation within a photosynthetic eukaryote.

Epithemia spp. diatoms contain obligate, nitrogen-fixing endosymbionts, or "diazoplasts", derived from cyanobacteria. These algae are a rare example of photosynthetic eukaryotes that have successfully coupled oxygenic photosynthesis with oxygen-sensitive nitrogenase activity. Here, we report a newly-isolated species, E. clementina, as a model to investigate endosymbiotic acquisition of nitrogen fixation. To detect the metabolic changes associated with endosymbiotic specialization, we compared nitrogen fixation, associated carbon and nitrogen metabolism, and their regulatory pathways in the Epithemia diazoplast with its close, free-living cyanobacterial relative, Crocosphaera subtropica. Unlike C. subtropica, we show that nitrogenase activity in the diazoplast is concurrent with, and even dependent on, host photosynthesis and no longer associated with cyanobacterial glycogen storage suggesting carbohydrates are imported from the host diatom. Carbohydrate catabolism in the diazoplast indicates that the oxidative pentose pathway and oxidative phosphorylation, in concert, generates reducing equivalents and ATP and consumes oxygen to support nitrogenase activity. In contrast to expanded nitrogenase activity, the diazoplast has diminished ability to utilize alternative nitrogen sources. Upon ammonium repletion, negative feedback regulation of nitrogen fixation was conserved, however ammonia assimilation showed paradoxical responses in the diazoplast compared with C. subtropica. The altered nitrogen regulation likely favors nitrogen transfer to the host. Our results suggest that the diazoplast is specialized for endosymbiotic nitrogen fixation. Altogether, we establish a new model for studying endosymbiosis, perform the first functional characterization of this diazotroph endosymbiosis, and identify metabolic adaptations for endosymbiotic acquisition of a critical biological function.

Comparing the transmission potential from sequence and surveillance data of 2009 North American influenza pandemic waves.

Technological advancements in phylodynamic modeling coupled with the accessibility of real-time pathogen genetic data are increasingly important for understanding the infectious disease transmission dynamics. In this study, we compare the transmission potentials of North American influenza A(H1N1)pdm09 derived from sequence data to that derived from surveillance data. The impact of the choice of tree-priors, informative epidemiological priors, and evolutionary parameters on the transmission potential estimation is evaluated. North American Influenza A(H1N1)pdm09 hemagglutinin (HA) gene sequences are analyzed using the coalescent and birth-death tree prior models to estimate the basic reproduction number (R 0 ). Epidemiological priors gathered from published literature are used to simulate the birth-death skyline models. Path-sampling marginal likelihood estimation is conducted to assess model fit. A bibliographic search to gather surveillance-based R 0 values were consistently lower (mean ≤ 1.2) when estimated by coalescent models than by the birth-death models with informative priors on the duration of infectiousness (mean ≥ 1.3 to ≤2.88 days). The user-defined informative priors for use in the birth-death model shift the directionality of epidemiological and evolutionary parameters compared to non-informative estimates. While there was no certain impact of clock rate and tree height on the R 0 estimation, an opposite relationship was observed between coalescent and birth-death tree priors. There was no significant difference (p = 0.46) between the birth-death model and surveillance R 0 estimates. This study concludes that tree-prior methodological differences may have a substantial impact on the transmission potential estimation as well as the evolutionary parameters. The study also reports a consensus between the sequence-based R 0 estimation and surveillance-based R 0 estimates. Altogether, these outcomes shed light on the potential role of phylodynamic modeling to augment existing surveillance and epidemiological activities to better assess and respond to emerging infectious diseases.

Direct long read visualization reveals metabolic interplay between two antimalarial drug targets.

Increases in the copy number of large genomic regions, termed genome amplification, are an important adaptive strategy for malaria parasites. Numerous amplifications across the Plasmodium falciparum genome contribute directly to drug resistance or impact the fitness of this protozoan parasite. During the characterization of parasite lines with amplifications of the dihydroorotate dehydrogenase (DHODH) gene, we detected increased copies of an additional genomic region that encompassed 3 genes (~5 kb) including GTP cyclohydrolase I (GCH1 amplicon). While this gene is reported to increase the fitness of antifolate resistant parasites, GCH1 amplicons had not previously been implicated in any other antimalarial resistance context. Here, we further explored the association between GCH1 and DHODH copy number. Using long read sequencing and single read visualization, we directly observed a higher number of tandem GCH1 amplicons in parasites with increased DHODH copies (up to 9 amplicons) compared to parental parasites (3 amplicons). While all GCH1 amplicons shared a consistent structure, expansions arose in 2-unit steps (from 3 to 5 to 7, etc copies). Adaptive evolution of DHODH and GCH1 loci was further bolstered when we evaluated prior selection experiments; DHODH amplification was only successful in parasite lines with pre-existing GCH1 amplicons. These observations, combined with the direct connection between metabolic pathways that contain these enzymes, lead us to propose that the GCH1 locus is beneficial for the fitness of parasites exposed to DHODH inhibitors. This finding highlights the importance of studying variation within individual parasite genomes as well as biochemical connections of drug targets as novel antimalarials move towards clinical approval.

Mitochondrial genomic data are effective at placing mycoheterotrophic lineages in plant phylogeny.

Fully mycoheterotrophic plants can be difficult to place in plant phylogeny due to elevated substitution rates associated with photosynthesis loss. This potentially limits the effectiveness of downstream analyses of mycoheterotrophy that depend on accurate phylogenetic inference. Although mitochondrial genomic data sets are rarely used in plant phylogenetics, theory predicts that they should be resilient to long-branch artefacts, thanks to their generally slow evolution, coupled with limited rate elevation in heterotrophs. We examined the utility of mitochondrial genomes for resolving contentious higher-order placements of mycoheterotrophic lineages in two test cases: monocots (focusing on Dioscoreales) and Ericaceae. We find Thismiaceae to be distantly related to Burmanniaceae in the monocot order Dioscoreales, conflicting with current classification schemes based on few gene data sets. We confirm that the unusual Afrothismia is related to Taccaceae-Thismiaceae, with a corresponding independent loss of photosynthesis. In Ericaceae we recovered the first well supported relationships among its five major lineages: mycoheterotrophic Ericaceae are not monophyletic, as pyroloids are inferred to be sister to core Ericaceae, and monotropoids to arbutoids. Genes recovered from mitochondrial genomes collectively resolved previously ambiguous mycoheterotroph higher-order relationships. We propose that mitochondrial genomic data should be considered in standardised gene panels for inferring overall plant phylogeny.

Assessment of stream macroinvertebrate communities with eDNA is not congruent with tissue-based metabarcoding.

Freshwater biomonitoring programmes routinely sample aquatic macroinvertebrates. These samples are time-consuming to collect, as well as challenging and costly to identify reliably genus or species. Environmental DNA (eDNA) metabarcoding has emerged as a surrogate to traditional collection techniques and has been used in whole-community approaches across several taxa and ecosystems. However, the usefulness of eDNA-based detection of freshwater macroinvertebrates has not been extensively explored. Few studies have directly compared bulk sample and eDNA metabarcoding at a local scale to assess how effective each method is at characterizing aquatic macroinvertebrate communities. Here, we collected both eDNA and kicknet samples at the same sample transect locations across nine different streams in southern Ontario, Canada. We observed minimal overlap in community composition between these paired samples. Bulk tissue metabarcoding resulted in a greater proportion of sequences belonging to metazoan taxa (over 99%) than eDNA (12%) and had higher OTU richness for macroinvertebrate taxa. We suggest that degenerate primers are not effective for eDNA metabarcoding due to the high degree of nontarget amplification and subsequently low yield of target DNA. While both bulk sample and eDNA metabarcoding had the power to detect differences between stream communities, eDNA did not represent local communities. Bulk tissue metabarcoding thus provides a more accurate representation of local stream macroinvertebrate communities and is the preferred method if smaller-scale spatial resolution is an important factor in data analyses.

Validation of COI metabarcoding primers for terrestrial arthropods.

Metabarcoding can rapidly determine the species composition of bulk samples and thus aids biodiversity and ecosystem assessment. However, it is essential to use primer sets that minimize amplification bias among taxa to maximize species recovery. Despite this fact, the performance of primer sets employed for metabarcoding terrestrial arthropods has not been sufficiently evaluated. This study tests the performance of 36 primer sets on a mock community containing 374 insect species. Amplification success was assessed with gradient PCRs and the 21 most promising primer sets selected for metabarcoding. These 21 primer sets were also tested by metabarcoding a Malaise trap sample. We identified eight primer sets, mainly those including inosine and/or high degeneracy, that recovered more than 95% of the species in the mock community. Results from the Malaise trap sample were congruent with the mock community, but primer sets generating short amplicons produced potential false positives. Taxon recovery from both mock community and Malaise trap sample metabarcoding were used to select four primer sets for additional evaluation at different annealing temperatures (40-60 °C) using the mock community. The effect of temperature varied by primer pair but overall it only had a minor effect on taxon recovery. This study reveals the weak performance of some primer sets employed in past studies. It also demonstrates that certain primer sets can recover most taxa in a diverse species assemblage. Thus, based our experimental set up, there is no need to employ several primer sets targeting the same gene region. We identify several suitable primer sets for arthropod metabarcoding, and specifically recommend BF3 + BR2, as it is not affected by primer slippage and provides maximal taxonomic resolution. The fwhF2 + fwhR2n primer set amplifies a shorter fragment and is therefore ideal when targeting degraded DNA (e.g., from gut contents).

Metabarcoding a diverse arthropod mock community.

Although DNA metabarcoding is an attractive approach for monitoring biodiversity, it is often difficult to detect all the species present in a bulk sample. In particular, sequence recovery for a given species depends on its biomass and mitome copy number as well as the primer set employed for PCR. To examine these variables, we constructed a mock community of terrestrial arthropods comprised of 374 species. We used this community to examine how species recovery was impacted when amplicon pools were constructed in four ways. The first two protocols involved the construction of bulk DNA extracts from different body segments (Bulk Abdomen, Bulk Leg). The other protocols involved the production of DNA extracts from single legs which were then merged prior to PCR (Composite Leg) or PCR-amplified separately (Single Leg) and then pooled. The amplicons generated by these four treatments were then sequenced on three platforms (Illumina MiSeq, Ion Torrent PGM and Ion Torrent S5). The choice of sequencing platform did not substantially influence species recovery, although the Miseq delivered the highest sequence quality. As expected, species recovery was most efficient from the Single Leg treatment because amplicon abundance varied little among taxa. Among the three treatments where PCR occurred after pooling, the Bulk Abdomen treatment produced a more uniform read abundance than the Bulk Leg or Composite Leg treatment. Primer choice also influenced species recovery and evenness. Our results reveal how variation in protocols can have substantial impacts on perceived diversity unless sequencing coverage is sufficient to reach an asymptote.

Convergent Plastome Evolution and Gene Loss in Holoparasitic Lennoaceae.

The Lennoaceae, a small monophyletic plant family of root parasites endemic to the Americas, are one of the last remaining independently evolved lineages of parasitic angiosperms lacking a published plastome. In this study, we present the assembled and annotated plastomes of two species spanning the crown node of Lennoaceae, Lennoa madreporoides and Pholisma arenarium, as well as their close autotrophic relative from the sister family Ehretiaceae, Tiquilia plicata. We find that the plastomes of L. madreporoides and P. arenarium are similar in size and gene content, and substantially reduced compared to T. plicata, consistent with trends seen in other holoparasitic lineages. In particular, most plastid genes involved in photosynthesis function have been lost, whereas housekeeping genes (ribosomal protein-coding genes, rRNAs, and tRNAs) are retained. One notable exception is the persistence of a rbcL open reading frame in P. arenarium but not L. madreporoides suggesting a nonphotosynthetic function for this gene. Of the retained coding genes, dN/dS ratios indicate that some remain under purifying selection, whereas others show relaxed selection. Overall, this study supports the mounting evidence for convergent plastome evolution in flowering plants following the shift to heterotrophy.

Finding the pond through the weeds: eDNA reveals underestimated diversity of pondweeds.

PREMISE OF THE STUDY: The detection of environmental DNA (eDNA) using high-throughput sequencing has rapidly emerged as a method to detect organisms from environmental samples. However, eDNA studies of aquatic biomes have focused on surveillance of animal species with less emphasis on plants. Pondweeds are important bioindicators of freshwater ecosystems, although their diversity is underestimated due to difficulties in morphological identification and monitoring. METHODS: A protocol was developed to detect pondweeds in water samples using atpB-rbcL and ITS2 markers. The water samples were collected from the Grand River within the rare Charitable Research Reserve, Ontario (RARE). Short fragments were amplified using primers targeting pondweeds, sequenced on an Ion Torrent Personal Genome Machine, and assigned to the taxonomy using a local DNA reference library and GenBank. RESULTS: We detected two species earlier documented at the experimental site during ecological surveys (Potamogeton crispus and Stuckenia pectinata) and three species new to the RARE checklist (P. foliosus, S. filiformis, and Zannichellia palustris). DISCUSSION: Our targeted approach to track the species composition of pondweeds in freshwater ecosystems revealed underestimation of their diversity. This result suggests that eDNA is an effective tool for monitoring plant diversity in aquatic habitats.

A Sequel to Sanger: amplicon sequencing that scales.

BACKGROUND: Although high-throughput sequencers (HTS) have largely displaced their Sanger counterparts, the short read lengths and high error rates of most platforms constrain their utility for amplicon sequencing. The present study tests the capacity of single molecule, real-time (SMRT) sequencing implemented on the SEQUEL platform to overcome these limitations, employing 658 bp amplicons of the mitochondrial cytochrome c oxidase I gene as a model system. RESULTS: By examining templates from more than 5000 species and 20,000 specimens, the performance of SMRT sequencing was tested with amplicons showing wide variation in GC composition and varied sequence attributes. SMRT and Sanger sequences were very similar, but SMRT sequencing provided more complete coverage, especially for amplicons with homopolymer tracts. Because it can characterize amplicon pools from 10,000 DNA extracts in a single run, the SEQUEL can reduce greatly reduce sequencing costs in comparison to first (Sanger) and second generation platforms (Illumina, Ion). CONCLUSIONS: SMRT analysis generates high-fidelity sequences from amplicons with varying GC content and is resilient to homopolymer tracts. Analytical costs are low, substantially less than those for first or second generation sequencers. When implemented on the SEQUEL platform, SMRT analysis enables massive amplicon characterization because each instrument can recover sequences from more than 5 million DNA extracts a year.

Escaping introns in COI through cDNA barcoding of mushrooms: Pleurotus as a test case.

DNA barcoding involves the use of one or more short, standardized DNA fragments for the rapid identification of species. A 648-bp segment near the 5' terminus of the mitochondrial cytochrome c oxidase subunit I (COI) gene has been adopted as the universal DNA barcode for members of the animal kingdom, but its utility in mushrooms is complicated by the frequent occurrence of large introns. As a consequence, ITS has been adopted as the standard DNA barcode marker for mushrooms despite several shortcomings. This study employed newly designed primers coupled with cDNA analysis to examine COI sequence diversity in six species of Pleurotus and compared these results with those for ITS. The ability of the COI gene to discriminate six species of Pleurotus, the commonly cultivated oyster mushroom, was examined by analysis of cDNA. The amplification success, sequence variation within and among species, and the ability to design effective primers was tested. We compared ITS sequences to their COI cDNA counterparts for all isolates. ITS discriminated between all six species, but some sequence results were uninterpretable, because of length variation among ITS copies. By comparison, a complete COI sequences were recovered from all but three individuals of Pleurotus giganteus where only the 5' region was obtained. The COI sequences permitted the resolution of all species when partial data was excluded for P. giganteus. Our results suggest that COI can be a useful barcode marker for mushrooms when cDNA analysis is adopted, permitting identifications in cases where ITS cannot be recovered or where it offers higher resolution when fresh tissue is. The suitability of this approach remains to be confirmed for other mushrooms.

On the brink: the highly reduced plastomes of nonphotosynthetic Ericaceae.

Ericaceae (the heather family) is a large and diverse group of plants that forms elaborate symbiotic relationships with mycorrhizal fungi, and includes several nonphotosynthetic lineages. Using an extensive sample of fully mycoheterotrophic (MH) species, we explored inter- and intraspecific variation as well as selective constraints acting on the plastomes of these unusual plants. The plastomes of seven MH genera were analysed in a phylogenetic context with two geographically disparate individuals sequenced for Allotropa, Monotropa, and Pityopus. The plastomes of nonphotosynthetic Ericaceae are highly reduced in size (c. 33-41 kbp) and content, having lost all photosynthesis-related genes, and are reduced to encoding housekeeping genes as well as a protease subunit (clpP)-like and acetyl-CoA carboxylase subunit D (accD)-like open reading frames. Despite an increase in the rate of their nucleotide substitutions, the remaining protein-coding genes are typically under purifying selection in full MHs. We also identified ribosomal proteins under relaxed or neutral selection. These plastomes also exhibit striking structural rearrangements. Intraspecific variation within MH Ericaceae ranges from a few differences (Allotropa) to extensive population divergences (Monotropa, Hypopitys), which indicates that cryptic speciation may be occurring in several lineages. The pattern of gene loss within fully MH Ericaceae plastomes suggests an advanced state of degradation.

Testing the Efficacy of DNA Barcodes for Identifying the Vascular Plants of Canada.

Their relatively slow rates of molecular evolution, as well as frequent exposure to hybridization and introgression, often make it difficult to discriminate species of vascular plants with the standard barcode markers (rbcL, matK, ITS2). Previous studies have examined these constraints in narrow geographic or taxonomic contexts, but the present investigation expands analysis to consider the performance of these gene regions in discriminating the species in local floras at sites across Canada. To test identification success, we employed a DNA barcode reference library with sequence records for 96% of the 5108 vascular plant species known from Canada, but coverage varied from 94% for rbcL to 60% for ITS2 and 39% for matK. Using plant lists from 27 national parks and one scientific reserve, we tested the efficacy of DNA barcodes in identifying the plants in simulated species assemblages from six biogeographic regions of Canada using BLAST and mothur. Mean pairwise distance (MPD) and mean nearest taxon distance (MNTD) were strong predictors of barcode performance for different plant families and genera, and both metrics supported ITS2 as possessing the highest genetic diversity. All three genes performed strongly in assigning the taxa present in local floras to the correct genus with values ranging from 91% for rbcL to 97% for ITS2 and 98% for matK. However, matK delivered the highest species discrimination (~81%) followed by ITS2 (~72%) and rbcL (~44%). Despite the low number of plant taxa in the Canadian Arctic, DNA barcodes had the least success in discriminating species from this biogeographic region with resolution ranging from 36% with rbcL to 69% with matK. Species resolution was higher in the other settings, peaking in the Woodland region at 52% for rbcL and 87% for matK. Our results indicate that DNA barcoding is very effective in identifying Canadian plants to a genus, and that it performs well in discriminating species in regions where floristic diversity is highest.

Using herbarium-derived DNAs to assemble a large-scale DNA barcode library for the vascular plants of Canada.

PREMISE OF THE STUDY: Constructing complete, accurate plant DNA barcode reference libraries can be logistically challenging for large-scale floras. Here we demonstrate the promise and challenges of using herbarium collections for building a DNA barcode reference library for the vascular plant flora of Canada. METHODS: Our study examined 20,816 specimens representing 5076 of 5190 vascular plant species in Canada (98%). For 98% of the specimens, at least one of the DNA barcode regions was recovered from the plastid loci rbcL and matK and from the nuclear ITS2 region. We used beta regression to quantify the effects of age, type of preservation, and taxonomic affiliation (family) on DNA sequence recovery. RESULTS: Specimen age and method of preservation had significant effects on sequence recovery for all markers, but influenced some families more (e.g., Boraginaceae) than others (e.g., Asteraceae). DISCUSSION: Our DNA barcode library represents an unparalleled resource for metagenomic and ecological genetic research working on temperate and arctic biomes. An observed decline in sequence recovery with specimen age may be associated with poor primer matches, intragenomic variation (for ITS2), or inhibitory secondary compounds in some taxa.

Correction: Authentication of Herbal Supplements Using Next-Generation Sequencing.

[This corrects the article DOI: 10.1371/journal.pone.0156426.].

Authentication of Herbal Supplements Using Next-Generation Sequencing.

BACKGROUND: DNA-based testing has been gaining acceptance as a tool for authentication of a wide range of food products; however, its applicability for testing of herbal supplements remains contentious. METHODS: We utilized Sanger and Next-Generation Sequencing (NGS) for taxonomic authentication of fifteen herbal supplements representing three different producers from five medicinal plants: Echinacea purpurea, Valeriana officinalis, Ginkgo biloba, Hypericum perforatum and Trigonella foenum-graecum. Experimental design included three modifications of DNA extraction, two lysate dilutions, Internal Amplification Control, and multiple negative controls to exclude background contamination. Ginkgo supplements were also analyzed using HPLC-MS for the presence of active medicinal components. RESULTS: All supplements yielded DNA from multiple species, rendering Sanger sequencing results for rbcL and ITS2 regions either uninterpretable or non-reproducible between the experimental replicates. Overall, DNA from the manufacturer-listed medicinal plants was successfully detected in seven out of eight dry herb form supplements; however, low or poor DNA recovery due to degradation was observed in most plant extracts (none detected by Sanger; three out of seven-by NGS). NGS also revealed a diverse community of fungi, known to be associated with live plant material and/or the fermentation process used in the production of plant extracts. HPLC-MS testing demonstrated that Ginkgo supplements with degraded DNA contained ten key medicinal components. CONCLUSION: Quality control of herbal supplements should utilize a synergetic approach targeting both DNA and bioactive components, especially for standardized extracts with degraded DNA. The NGS workflow developed in this study enables reliable detection of plant and fungal DNA and can be utilized by manufacturers for quality assurance of raw plant materials, contamination control during the production process, and the final product. Interpretation of results should involve an interdisciplinary approach taking into account the processes involved in production of herbal supplements, as well as biocomplexity of plant-plant and plant-fungal biological interactions.

Plastid genome evolution across the genus Cuscuta (Convolvulaceae): two clades within subgenus Grammica exhibit extensive gene loss.

The genus Cuscuta (Convolvulaceae, the morning glory family) is one of the most intensely studied lineages of parasitic plants. Whole plastome sequencing of four Cuscuta species has demonstrated changes to both plastid gene content and structure. The presence of photosynthetic genes under purifying selection indicates that Cuscuta is cryptically photosynthetic. However, the tempo and mode of plastid genome evolution across the diversity of this group (~200 species) remain largely unknown. A comparative investigation of plastid genome content, grounded within a phylogenetic framework, was conducted using a slot-blot Southern hybridization approach. Cuscuta was extensively sampled (~56% of species), including groups previously suggested to possess more altered plastomes compared with other members of this genus. A total of 56 probes derived from all categories of protein-coding genes, typically found within the plastomes of flowering plants, were used. The results indicate that two clades within subgenus Grammica (clades 'O' and 'K') exhibit substantially more plastid gene loss relative to other members of Cuscuta. All surveyed members of the 'O' clade show extensive losses of plastid genes from every category of genes typically found in the plastome, including otherwise highly conserved small and large ribosomal subunits. The extent of plastid gene losses within this clade is similar in magnitude to that observed previously in some non-asterid holoparasites, in which the very presence of a plastome has been questioned. The 'K' clade also exhibits considerable loss of plastid genes. Unlike in the 'O' clade, in which all species seem to be affected, the losses in clade 'K' progress phylogenetically, following a pattern consistent with the Evolutionary Transition Series hypothesis. This clade presents an ideal opportunity to study the reduction of the plastome of parasites 'in action'. The widespread plastid gene loss in these two clades is hypothesized to be a consequence of the complete loss of photosynthesis. Additionally, taxa that would be the best candidates for entire plastome sequencing are identified in order to investigate further the loss of photosynthesis and reduction of the plastome within Cuscuta.

Plastid genome evolution in mycoheterotrophic Ericaceae.

Unlike parasitic plants, which are linked to their hosts directly through haustoria, mycoheterotrophic (MHT) plants derive all or part of their water and nutrients from autothrophs via fungal mycorrhizal intermediaries. Ericaceae, the heather family, are a large and diverse group of plants known to form elaborate symbiotic relationships with mycorrhizal fungi. Using PHYA sequence data, we first investigated relationships among mycoheterotrophic Ericaceae and their close autotrophic relatives. Phylogenetic results suggest a minimum of two independent origins of MHT within this family. Additionally, a comparative investigation of plastid genomes (plastomes) grounded within this phylogenetic framework was conducted using a slot-blot Southern hybridization approach. This survey encompassed numerous lineages of Ericaceae with different life histories and trophic levels, including multiple representatives from mixotrophic Pyroleae and fully heterotrophic Monotropeae and Pterosporeae. Fifty-four probes derived from all categories of protein coding genes typically found within the plastomes of flowering plants were used. Our results indicate that the holo-mycoheterotrophic Ericaceae exhibit extensive loss of genes relating to photosynthetic function and expression of the plastome but retain genes with possible functions outside photosynthesis. Mixotrophic taxa tend to retain most genes relating to photosynthetic functions but are varied regarding the plastid ndh gene content. This investigation extends previous inferences that the loss of the NDH complex occurs prior to becoming holo-heterotrophic and it shows that the pattern of gene losses among mycoheterotrophic Ericaceae is similar to that of haustorial parasites. Additionally, we identify the most desirable candidate species for entire plastome sequencing.

Loss of all plastid ndh genes in Gnetales and conifers: extent and evolutionary significance for the seed plant phylogeny.

The exact phylogenetic position of Gnetales, a small, highly modified group of gymnosperms with an accelerated rate of molecular evolution, is one of the most challenging issues for seed plant systematics. Recent results from entire plastid genome (ptDNA) sequencing revealed the absence of the entire suite of plastid ndh genes in several species of Gnetales and the pine family (Pinaceae) potentially highlighting a major structural feature linking these two groups-concerted loss of all plastid genes for the NADH dehydrogenase complex. However, the precise extent of ndh gene loss in gymnosperms has not been surveyed. Using a slot-blot hybridization method, we probed all 11 ndh genes in 162 species from 70 of 85 gymnosperm genera. We find that all ndh genes are absent across Gnetales and Pinaceae, but not in any other group of gymnosperms. This feature represents either a major synapomorphy for a clade consisting of these two lineages or, less likely, a convergent loss. Our survey substantially extends previous inferences based on more limited sampling and, if the former evolutionary interpretation is correct, it provides additional support for the contentious "gnepine" hypothesis, which places Gnetales as sister to Pinaceae.